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Abstract: 

Coupled double well {(p"^) one-dimensional potentials abound in both condensed matter physics and 
field theory. Here we provide an exhaustive set of exact periodic solutions of a coupled (j)^ model in an 
external field in terms of elliptic functions (domain wall arrays) and obtain single domain wall solutions in 
specific limits. We also calculate the energy and interaction between solitons for various solutions. Both 
topological and nontopological (e.g. some pulse-like solutions in the presence of a conjugate field) domain 
walls are obtained. We relate some of these solutions to the recently observed magnetic domain walls in 
certain multiferroic materials and also in the field theory context wherever possible. Discrete analogs of 
these coupled models, relevant for structural transitions on a lattice, are also considered. 
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1 Introduction 



There are many physical situations, both in condensed matter and field theory, where two double well 
potentials model the phenomena of interest with a specific coupling between the two fields. One such 
phenomenon of current intense interest is the coexistence of magnetism and ferroelectricity (i.e. magneto- 
electricity) in a given material. This is a highly desired functionality in technological applications involving 
cross-field response, switching and actuation. In general, this phenomenon is referred to as multiferroic 
behavior PP. Recently, two different classes of (single phase) multiferroics, namely the orthorhombically 
distorted perovskites ^ and rare earth hexagonal structures |3j, have emerged. The latter show magnetic 
domain walls in the basal planes which can be modeled by a coupled (j)^ model ^ in the presence of a 
magnetic field. Coupled (p^ models [SUHIIZ] also arise in the context of many ferroelectric and other second 
order phase transitions. The coupled (p^ model for multiferroics jlj has a biquadratic coupling whereas the 
coupled (/)^ model for a surface phase transition with hydration forces [7j, relevant in biophysics context, 
has a bilinear coupling. Other types of couplings are also known for structural phase transitions with strain 

Similarly, there are analogous coupled models in field theory [^1^]. Several related models have been 
discussed in the literature and their soliton solutions have been found ^1 El El El El El including 
periodic ones |191 12UI I21j . Here our motivation is to obtain various possible domain wall solutions of these 
models with either a bilinear or a biquadratic coupling and then connect to experimental observations 
wherever possible. 

The paper is organized as follows. In Sec. II we provide the solutions for the coupled <j)^ model with 
an explicit biquadratic coupling in the presence of an external field (with an additional linear-quadratic 
coupling) and calculate their energy as well as interaction between the solitons. We also obtain solutions 
in the limit of no field. In Sec. Ill we consider a coupled 0^ model with a bilinear coupling. In Sec. IV we 
obtain several solutions of a coupled discrete 0^ model with biquadratic coupling (but without an external 
field). In Sec. V we obtain a solution of the coupled discrete (/>^ model with bilinear coupling. To the best 
of our knowledge, no solution of the coupled t/)^ model in an external field is known. Similarly, no solutions 
of either the coupled continuum (p'^ model with a bilinear coupling or the discrete case are known. Even 
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for the continuum coupled case with a biquadratic coupUng and zero external field, out of the six possible 
solutions only three were known previously ^20, but the other three are new. Finally, we conclude in Sec. 
VI with remarks on related models. 



2 Coupled (j)^ solutions in an external field 

We consider several exact solutions of the coupled (f>^ system in a magnetic field (Hz) as given in Ref. 
[1] for hexagonal multiferroics. In particular, there are nine periodic solutions (valid for arbitrary m, the 
modulus of elliptic functions), which at m = 1 reduce to just four solutions. In particular, there is one 
"bright-bright", one "dark-dark" and one each of bright-dark and dark-bright solutions. Notice that the 
equations of motion are asymmetric in the two scalar fields {(j) and ip) due to different coupling of the scalar 
fields to the magnetic field. Thus, the dark-bright and the bright-dark solutions are distinct. 

The potential, with a biquadratic coupling between the two fields and in an external magnetic field 
(i?^), is given by 

where Oj, Pi, 7 and pi are material (or system) dependent parameters. Hence the (static) equations of 
motion are 

= 2aicP + 4/3i(/.3 + 270V' - H,[pi + 3p2(p^ + ps^P^] , (2) 
^ = 2a2^p + 4/32^'^ + 27</'V - 2H,p3cj)iP . (3) 

These coupled set of equations admit several periodic solutions which we now discuss one by one system- 
atically. 

For static solutions the energy is given by 
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E 



dx , (4) 



where the limits of integration are from —00 to 00 in the case of hyperbolic solutions (i.e., single solitons) 
on the full line. On the other hand, in the case of periodic solutions (i.e. soliton lattices), the limits are 
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from —K{m) to +K{m). Here K{m) [and E{m) below] denote the complete elliptic integral of the first 
(and second) kind Using equations of motion, one can show that for all of our solutions 



1 (<^\\\(d±\ 



2 \dx J ~^ 2 \dx J 



+ C. 



(5) 



where the constant C in general varies from solution to solution. Hence the energy E = E — J C dx is 
given by 



E = E- 



Cdx 



dx. 



(6) 



^dx J \dx , 

Below we will give explicit expressions for energy in the case of all nine periodic solutions (and hence the 
corresponding four hyperbolic solutions). In each case we also provide an expression for the constant C. 

2.1 Solution I 

We look for the most general periodic solutions in terms of the Jacobi elliptic functions sn(x, m), cn{x, m) 
and dn(x,m) It is not difficult to show that 



4> = F + Asn[D{x + xo),m] , ip = G + Bsn[D{x + xo),m] , 

is an exact solution provided the following eight coupled equations are satisfied 

2aiF + 4/3iF3 + 2jFG^ - H^pi - 3H,p2F^ - H.psG^ = , 

2aiA + l2(3iF'^A + A-yBFG + 2-fAG'^ - 6H^p2AF - 2H^p^BG = -(1 + m)AD'^ 
12(3iFA^ + 2-fFB'^ + A-fABG - 3H^p2A'^ - if^psB^ = , 
2PiA^ + = mD^ , 
202^ + 4/32^3 + 2-fGF^ - 2H,p3GF = , 
202^ + 12/52^^5 + A'yAFG + 2-fBF^ - 2H^p^{BF + AG) = -(1 + m)BD'^ , 
12I32GB^ + 2-iGA^ + A^iABF - 2H^p^AB = , 



(7) 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 
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Here A and B denote the amplitudes of the "kink lattice", F and G are constants, D is an inverse 
characteristic length and xq is the (arbitrary) location of the kink. Five of these equations determine 
the five unknowns A,B,D,F,G while the other three equations, give three constraints between the nine 
parameters ai^2, (!^i,2,'y, H^, pi, P2, Ps- In particular, A and B are given by 

,_ mD^[2(32-7] 2_ mD^Wi-l] .... 
[4/3i/?2-7'] ' [4/31/32 -7^] ' ^ ' 

Thus this solution exists provided 2/?i > 7 and 2/?2 > 7. 

It may be noted here that in case both F,G = ^ then no solution exists so long as Hz 7^ 0. In fact this 

is true for all the nine solutions that we discuss below. There are two special cases (i.e., when either F = 

or G = 0) when the analysis becomes somewhat simpler. 

G=0, Fy^O: 

In this case A, B are again given by Eq. (|16l) while 

(l + m)Z?2 = ^^-^, F=^. (17) 
27 p-i 27 

The three (and not four, since one of the equations is identically satisfied) constraints are 

m P2P3HI pa pjHl P2P3HI pi'y 

Pi = 7^ , ai = — ^ H > "2 = -. + • (18) 

2/03 27 ps 47 47 p3 

m=l: 

In this limiting case we have a bright-bright soliton solution given by 

(/> = F + ^tanh[L'(x + xo)] , -0 = 5tanh[L'(x + xq)] , (19) 
with A, B and D determined by 

2_ D^m-i] 2_ D^m-i] .p2 _ P2P3H! 2p,j 
mh - 7^] ' mh - 7^] ' 27 P3 ' ^ ^ 

while the other relations remain unchanged and are again given by Eqs. H17() and ()18() . 
F=0, G/0: 

In this case A, B are again given by Eq. H16|) with 

(l + n.)i^^ = 4a2 + ^^^^, ^' = -^' 
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and the corresponding four constraints are 

132 = — — , (1 + m)!)^ = -2ai - 27G^ H , 

2pi A 

i^ABG = m,p2A^ + H.psB^ , P3H,AB = -fOA^ + G/JsGS^ • (22) 

The solution at m = 1 can be easily written down as above. 
Special case of Hz = 0: 

In this case the field equations are completely symmetric in the two fields ^ and solution exists even 
if both F,G = In particular, with Hz = 0, the solution is as given by Eq. (7) but with F = Q = G [20] 
where A,B are again as given by Eq. (16) and further both ai,a2 turn out to be negative, i.e. 

{l + m)D^ 

01 = 02 = ^ . (23) 

Special case of 7^ = 4/3i/?2 

One can show that the solution ((71) exists even in case 7^ = 4/3i/32. It turns out such a solution exists 
only if 

2/?i = 2/32 = 7, (24) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

+ B^ = . (25) 
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Other relations can be easily worked out depending on if F or G (or neither) is zero. For example, in case 
G = 0, F / 0, one has 

P2 = P3, "2 = — , ai-a2 = ^^. (26) 

P3 27 

Energy: Corresponding to the periodic solution [Eq. Q with G = 0] the energy E and the constant C 
are given by 

E = ^i^^±^^ui ^ ^\E(m) - (1 - m)K(m)] , 
3m 

C = -^iA^ + B'^)D'^ - F\ai + WiF"^ - 2HzP2F] . (27) 

It is worth pointing out that even in the case of the solution ((TJ with either = 0, G 7^ or both F, G 
nonzero, the energy E is the same. Only the value of G is different. For example, in the case of F = 0, 



G 7^ 0, C is given by 



C = -\{A' + B^)D^-^, 



while in the case of both F, G being nonzero, C is 



On using the expansion formulas for E{m) and K{m) around m = 1 as given in 



(28) 



(29) 



K{m) = In 



+ 



(1 — m) 



E{m) = 1 + 



^Jl - mj ' 4 
(1 — m) 



In 



m 



+ ... 



In 



+ ..., 



(30) 



(31) 



."v/l — mJ 2 

for m near one, the energy of the periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (bright-bright) soliton solution [Eq. ()19() ] plus the interaction energy. We find 



E — Ekink + Eint — (^^ + B'^)D 



4 (1 — m) 

- + 

3 3 



(32) 



Note that this solution exists only when 2/3i > 7, 2/?2 > 7 and 4/3i/32 > 7^- The interaction energy vanishes 
at exactly m = 1, as it should. 



2.2 Solution II 



A different type of solution ("pulse lattice") is given by 



(f) = F + Acii[D{x + xo), , = G + Bcii[D{x + xq), "m-] , 
provided the following eight coupled equations are satisfied 

2aiF + A(3iF^ + 2^FG^ - H,pi - 3H,p2F^ - H^p^G^ = , 
2aiA + UPiF^A + 47^^^ + 27^^^ _ QH^p2AF - 2H^p3BG = (2m - 1)AD^ , 
UPiFA'^ + 2-fFB'^ + A-iABG - 3H,p2A'^ - H.psB"^ = , 
2piA^ + ^B^ = -mD^ , 
2a2G + 4p2G^ + 2-fGF^ - 2H,p3GF = , 



(33) 



(34) 

(35) 
(36) 
(37) 
(38) 



Thus this solution exists provided 7^ > 4/3i/3; 



202^ + 12I32G^B + A-yAFG + - 2H^p3{BF + AG) = (2m - 1)BD'^ , (39) 

12(32GB'^ + 2-/GA^ + A-/ABF - 2H^p-iAB = , (40) 

2/?2S^ + = -mD^ . (41) 

Notice that two of these equations are meaningful only if 7 < since /32 > from stability considerations. 
Thus, we write 7 = — 17|. Five of these equations determine the five unknowns A,B,D,F,G while other 
three equations, give three constraints between the nine parameters ai,2) /3i,2, 7, H^, pi, P2, P3- In particular, 
A and B are given by 

2 ^ mD'[2P2 + M 2 ^ mD^[2P, + \j\] 

[72-4/31/32] ' [7'-4/3i/32] ■ ^ ' 

•2- 

There are two special cases when the analysis becomes somewhat simpler and we consider both the 
cases one by one. 
G=0, Ft^O: 

In this case A, B are again given by Eq. 1)42(1 while 

{2m-l)D' = P^-^-^, F = -'-^, (43) 
^ ' 2|7| p3 2|7| ^ ' 

and the corresponding three constraints are 

2p3 ' ' 2|7| p, ' 4|7| + 4|7| p, ' 

F=0, G^O: 

In this case ^4, B are again given by Eq. (|42|) with 

(2m - 1)1)2 ^ - ^^^^ - 2|7|G2 , = -i| , (45) 
and the corresponding four constraints are 

a P3«2 -.^^2 . 2p3H^AG 

P2 = - — , 2m - 1)D = -4a2 5 , 

2pi B 

-4\-f\ABG = 3H,p2A^ + H,p3B\ p^H^AB = -\-f\GA^ + 6p2GB\ (46) 
Special case of Hz=0: 



In this case the field equations are completely symmetric in the two fields (/), ip and a solution exists 
even when both F,G = 0. In particular, with Hz = 0, the solution is as given by Eq. (|33jl but with 
F = = G where A,B are again as given by Eq. (|42j) and furthermore, ai,a2 are positive (negative) 
so long as m > (<) 1/2, i.e. 

ai = as = ^ . (47) 

Energy: Corresponding to the "pulse lattice" solution [Eq. (|c{3l) with G = 0] the energy is given by 

E = '—[{2m - l)E{m) + (1 - m)K{m)] , 

3m 

C = -^(1 - m){A^ + B^)D^ - F'^[ai + 3/?iF2 - 2HzP2F] . (48) 

It is worth pointing out that even in the case of the solution p3|) with either F = 0, G ^ or both F, 
G nonzero, the energy E is the same. Only the value of C is different. For example, in the case of F = 0, 
G / 0, C is given by 

C = -l{l-m){A' + B')D'-^, (49) 
while in the case of both F, G being nonzero, C is 

C = -^(1 - m)(^2 ^ ^2)^2 _ ^2[^^ ^ 3^^^2 _ 2HzP2F] - G^[a2 + S/JsG^ + 3jF^ - 2HzPzF] . (50) 

For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-dark) soliton solution 

(f) = F + Asec\i[D{x + xq)], = Bsec\i[D{x + xq)] , (51) 

plus the interaction energy. We find 

E = Epuise + Eint = (A'^ + B^)D 



2 5(1 -m) ,^ 

^ ^ + (1 - m) In 



(52) 



.3 6 V \/l — 

Note that this solution exists only when 7 < 0,4/?i/32 < 7^- Again, the interaction energy vanishes at 
m = 1. 
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2.3 Solution III 

In this case, there is another "pulse lattice" solution which is given by 

(j) = F + Adn[D{x + xo),m], ip = G + Bdn[D{x + xo),m] , (53) 
provided the following eight coupled equations are satisfied 

2aiF + 4/3iF3 + 2jFG^ - H,pi - 3H,p2F^ - H.psG'' = , (54) 

2aiA + l2(3iF^A + A'yBFG + 2jAG'^ - QH^p2AF - 2H^p3BG = (2 - m)AD'^ , (55) 

UPiFA^ + 2-fFB'^ + 4-fABG - m^p2A^ - H^p^B"^ = , (56) 

2l3iA^ + -iB'^ = -D^ , (57) 

2a2G + 4/?2G^ + 27GF2-2F^p3GF = 0, (58) 

2a2B + 12p2G'^B + A-fAFG + 2-fBF^ - 2H^p^{BF + AG) = (2 - m)BD'^ , (59) 

12^2GB^ + 2-iGA^ + 47A5F - 2H^p:iAB = , (60) 

2/32^2 + = -1)2 . (61) 

Notice that (as in the cn — cn case) two of these equations are meaningful only if 7 < since Pi,(32 > 
from stability considerations. We therefore write 7 = — 17|. Five of these equations determine the five 
unknowns A, B, D, F, G while the other three equations give three constraints between the nine parameters 
Q<i,2)/3i,2, h\, Hz, Pi, P2, P3- In particular, A and B are given by 

2 ^ D^m + i7i] „2 ^ D^m + M ..2^ 

[7^-4/31/32] ' [|7P-4/?i/?2]' ^ ^ 

There are two special cases when the analysis becomes somewhat simpler and we consider both the 
cases one by one. 
G=0, F/0: 

In this case A, B are again given by Eq. I|62|) while 

i2-m)D^ = P^-^-^, F = -^, (63) 
^ ' 2|7| p3 2|7| ^ ^ 
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and the three constraints are 

a P2\l\ ^ P2PZHI P1I7I ^ P2P3HI P1I7I 

PI = — 7, ' "1 = ' "2 = — -rr-r H tt-. ■ (D4j 

2/03 2|7| p3 4|7| 4|7| 

F=0, G/0: 

In this case A, B are again given by Eq. H62|) with 

(2-..)Z)^ = 2a,-^^^^-2|7|G^ G^ = -^, (65) 

and the corresponding four constraints are 

P2 = -7, — , 2 - m)D = -4a2 , 

2pi 5 

-4\^\ABG = 3H,p2A^ + H,p3B\ p^H^AB = -\^\GA^ + 6P2GB\ (66) 
Special case of Hz=0: 

With Hz = 0, the solution is as given by Eq. (|53|) but with F = = G ^20; where A, B are again as 
given by Eq. (|62() and furthermore, both ai,Q!2 are positive definite. 

(2-m)L»2 

ai = "2 = ^ . (67) 

Energy: Corresponding to the "pulse lattice" solution [Eq. with G = 0] the energy is given by 

E = ^ ^[(2 - m)E{m) - (1 - m)K{m)] , 

3 

C = ^(1 - m){A^ + ^2)1)2 _ F^^ai + 3/3iF2 - 2//,/)2F] . (68) 

It is worth pointing out that even in the case of the solution H53() with either F = 0, G ^ or both F, 
G nonzero, the energy E is the same. Only the value of C is different. For example, in the case of F = 0, 
G / 0, C is given by 

C = l(l-m)iA' + B')D'-^^, (69) 
while in the case of both F, G being nonzero, C is 

C = ^(1 - m)(^2 ^ ^2^^2 _ ^2[^^ ^ 3^^^2 _ 2HzP2F] - G'^[a2 + 3^20^ + 87^^ - 2HzP3F] . (70) 
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For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbohc (dark-dark) soliton solution as given by Eq. (|51jl plus the interaction energy. We find 

"2 (1-m) / 4 



E = Epulse + Eint = {A + B )D 



(1 — m) In ( , =\ 



(71) 



.3 2 

Note that this solution also exists only when 7 < 0, 4/3i/32 < 7^- Again, the interaction energy vanishes at 
m = 1. 

2.4 Solution IV 

In addition to the cn — cn and dn — dn solutions discussed above, there are two novel (mixed) soliton 
solutions of dn — cn and cn — dn type. Let us discuss them one by one. We shall see that for these two 
solutions (in fact it is true for all the six solutions that we discuss below) G is necessarily zero while F is 
necessarily nonzero (unless = 0), otherwise the solution does not exist. In particular, it is easily shown 
that 

<p = F + Adn[D{x + xo),m\ , = G + Bcn[D{x + xo),m\ , (72) 
is a solution provided G = and further, the following seven coupled equations are satisfied 

2aiF + 4piF^ + (2/m)7(m - 1)FB'^ - H,pi - SF^pa^^ + (1/m) (1 - m)H^p^B'^ = , (73) 

2aiA + UPiF'^A + (2/m)(m - 1)jAB'^ - 6H^p2AF = (2 - m)AD^ , (74) 

UPiFA^ + {2/m)jFB^ - 3H,p2A^ - {l/m)H,p^B'^ = , (75) 

7^2 + 2ml3iA^ = -niD"^ , (76) 

2a2B + 2(1 - m)-fA'^B + 27^^^ - 2H:,pj,BF = {2m - l)BD^ , (77) 

2p2B^ + m-/A^ = -mD^ , (78) 

A-fFAB - 2H^p3AB = . (79) 

Again a solution exists only if 7 < and thus we write 7 = — 17|. The solution (|72() exists provided 

2 ^ D'H+2P2] „2 ^ mD^[2/3i + |7|] ...^ 
[72-4/31/32] ' [7^ -4/31/32] ' ^ ^ 
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where 



(2 = 1(1 -,„)h|B=-?«M + ««!£i, F=-!!iEl, (81) 

m p3 2|7| 2|7| 



and the three constraints are 



o P2h\ P2P3HI P1I7I 

Pi = 7, ' "2 - 



(2m - 1)1)2 = 2q2- 2(1 -m)|7|^2^:^M. ^g2) 



2/33 2|7| P3 

■2I7I 

Special case of Hz = 0: 

With Hz = 0, the solution is as given by Eq. ((7^ but with F = = G where A and B are again as 
given by Eq. ()8U() and furthermore, 02 is positive, i.e. 

ai = — 2(1 - m)/?!^^ , ma2 = m— + 2(1 - m)(i2B'^ > . (83) 

Energy: Corresponding to the "pulse lattice" solution [Eq. ((7^ with G = 0] the energy E and C are 
given by 

2D / 

^ = — [(2 - m)mA^ + (2m - l)B'^]E{m) + (1 - m){B'^ - 2mA^)K{m] 
C = -A^D^{1 - m) - F^[ai + 3piF^ - 2HzP2F] + ^ ^[m-fF^ - agm + ^2(1 - m)B^]B^ .(84) 

For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbohc (dark-dark) soliton solution (|51|) plus the interaction energy. We find 



E — Epulse + Eint — D 



(85) 



.3 6 V \/l — m 

Note that this solution also exists only when 7 < 0, 4/3i/32 < 7^- The interaction energy vanishes for m = 1, 
as it should. 

2.5 Solution V 

It is easily shown that there is also a cn — dn solution which is distinct from the above dn — cn solution. 
This solution is given by 

(j) = F + Ac-D[D{x + XQ),m], V = G + -Bdn[L>(x + xq), m] , (86) 

13 



provided G = and the following seven coupled equations are satisfied 

2aiF + 4/?iF3 + 27(1 - m)FB^ - H^pi - 3H,p2F^ - (1 - m)H^p3B^ = , (87) 

2aiA + UpiF'^A + 2(1 - my/AB'^ - QH-,p2AF = -(1 - 2m)AD'^ , 

l2l3iFA^ + 2m-fFB^ - 3H^p2A^ - mH^paB^ = , 

mjB^ + 2piA^ = -mD^ , (90) 

2a2B - (2/m)(l - m)jA'^B + 2jBF'^ - 2H^p3BF = (2 - m)BD'^ , (91) 

2mp2B^ + 7^4^ = -mD^ , (92) 

A-fFAB - 2H^p3AB = . (93) 

Again a solution exists only if 7 < and hence we put 7 = — 17|. The solution turns out to be 

2 ^ mZ)2[|7|+2/?2] 2 ^ D'm + b\] .04^ 
[I7P-W2] ' [|7P-4/3i/32]' ^ ^ 



where 

2 

(9m. - ^\^)'^ = -9n - m\\^\R'^ , ^ _ 

P3 2|7| 2|7| 



/r, iNn2 on M I d2 2pi|7| P2P'iHl p^Hz 
(2m — l)D = — 2(1 — m)|7|i5 1 — — j — , F = — , (95) 



and the three constraints are 



a P2\l\ P2P?,Hl pl\ 

PI = — , ai ~ 



(2 - m)D2 = 2a2 + (2/m)(l - m)\-i\A^ + — . (96) 



2p3 2|7| p3 

2|7l 

Special case of Hz = 0: 

With i/^ = 0, the solution is as given by Eq. ()86() but with F = Q = G where A and B are again as 

given by Eq. ()94() and furthermore, ai is positive, i.e. 

r)2 mZ)2 
= m— + 2(1 - m)(3iA^ > , 02 = — 2(1 - 771)^2^^ ■ (97) 

It is worth pointing out that with = 0, the field equations are completely symmetric between the two 

fields (j) and ip and hence the solutions IV and V are identical in the limit = 0. 

m=l: 
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In the special case of m = 1 and G = 0, F ^ 0, all the four solutions II to V reduce to a dark-dark 
type soliton solution (|51|1 . i.e. 



F + Asech[D{x + xo)], = Bsech[D{x + xq)] , (98) 



with A, B, and D given by 



[I7P-W2]' [I7P-W2]' P3 2|7| ' ^ ^ 

while the other relations remain unchanged and are again given by Eqs. (|43() and (|44() . 

Similarly, in the special case of m = 1 and = 0, G 7^ 0, both the solutions II and III reduce to a 
dark-dark type soliton solution 

(/> = ^sech[L>(x + xo)] , ip = G + Bsech\D{x + xq)] , (100) 

with A and B again given by Eq. H99|) while 

D' = 2a^- '^P'^f^ -2\j\G\ = -4a2 - ^^2%^ . (101) 
A B 

Other relations remain unchanged and are again given by Eqs. H45|) and H46() . 

Energy: Corresponding to the "pulse lattice" solution [Eq. (|86|) with G = 0] the energy is given by 
2D / 

E = [(2 - m)mB'^ + (2m - l)^^]^;^^) + (1 - m){A'^ - 2mB'^)K{m)] 

3m V 

C = -^A^D\l - m) - F^[ai + SPiF^ - 2H,p2F] + [-7^^ + 03 + ^2(1 - m)B^]B^ . (102) 

For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-dark) soliton solution [Eq. (|51|) ] plus the interaction energy. We find 



E — Epuise + Eint — D 



2(^2 ^ ^2) ^ (L^(3ij2 _ 5^2) ^ _ ^)(^2 _ ^2) 1^ ^ 4 



(103) 



.3 6 — m, 

Note that this solution also exists only when 7 < 0, 4/3i/32 < 7^. Again, the interaction energy vanishes at 
m = 1 . It is amusing to note that the energy of the solution V is easily obtained from that of solution IV 
by simply interchanging A and B. 
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2.6 Solution VI 

Apart from the solutions which at m = 1 reduce to the bright-bright and dark-dark solutions, there are 

four solutions which at m = 1 go over to either bright-dark or dark-bright solutions, which we now discuss 
one by one. 

One such solution, kink-like in (j) and pulse-like in -0, is 

^ = F + A?,n[D{x + XQ),m], = G + Bcn[D{x + XQ),m] , (104) 

provided G = and the following seven coupled equations are satisfied 

2aiF + 4/3iF3 + 27^^^ - H^p^ - 3H,p2F^ - H.p^B^ = , (105) 

2aiA + uPiF'^A + IjAB"^ - 6H^p2AF = -(1 + m)AD'^ , (106) 

12/3iF^2 _ 2r^pQ2 _ sH,p2A^ + H.psB^ = , (107) 

2a2B + 2-fA^B + 2-fBF'^ - 2H^p^BF = (2m - 1)BD'^ , (109) 

2^2^^ - lA^ = -mD^ , (110) 

{'ijF-2H,p3)AB = 0. (Ill) 
In this case it turns out that the solution exists only if either 2/3i > 7 > 2/32, 7^ > 4/?i/?2 or 2/^2 > 7 > 



2/3i, Y < 4/3i/32- We find that 



2 _ mD'[^- 202] 2_ mD'm-^] 

[7'-4/3i/32] ' b^-mP2] ' ^ ' 



where 

2 _ o d2 , 2pi7 P2P3^^^ p _ PaHz 



-{l + m)D' = 2jB' + ^-^^^^, F = ^, (113) 

P3 ^7 ^7 



and the three constraints are 



P27 P2P3^z , Pil 

2p3 27 P3 

If 2 2 

{2m - 1)D^ = 2a2 + 27^12 _ _£i^3 _ ^1^4) 

27 
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Special case of = 0: 

With Hz = 0, the solution is again given by Eq. (|1U4|) but with F = Q = G where A and B are again 
as given by Eq. and furthermore, ai,a2 turn out to be negative, i.e. 

a, = J^±^-,B\a, = -^-2hB\ (115) 

Special case of 7^ = 4/3i/32 

One can show that the solution 1)104(1 exists even in case 7^^ = 4/3i/32- It turns out such a solution exists 
only if 

2/3i = 2/32=7, (116) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

-B^ = . (117) 

7 



Further, one has 



P2=Pz, ai = ^^ + ^, a2-ai = ^—-t\^. 118 
27 p3 2 27 



Energy: Corresponding to the mixed "kink-pulse lattice" solution [Eq. (|1U4|) ] the energy is given by 

E = — ([(2 - m)mB^ + (1 + m)A^]E{m) - (1 - m)iA^ + 2B'^)K(m) \ , 
3m V / 

C = -^A'^D^ - F^[ai + 3piF^ - 2H,p2F] + [-7^^ ^ ^ p^B^jB^ . (119) 

For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (bright-dark) soliton solution 

(j) = F + Atanh[D{x + xo)], = Bsech[D{x + xq)] , (120) 

plus the interaction energy. We find 



E — EsoHton + Eint — D 



:i2i) 



_3 6 V \/l — rn, 

Note that this solution also exists only when either 2(3i > 7 > 2(32 and 7^ > 4/3i/32 or 2/^2 > 7 > 2(3i and 
7^ < 4/3i/32- The interaction energy vanishes at m = 1. 
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2.7 Solution VII 

It is easy to show that another such (kink- and pulse-hke) solution is 

(t) = F + Asn[D{x + xo),m], = G + Bdn[D{x + xo),m\ , (122) 
provided G = and the foUowing seven coupled equations are satisfied 

2aiF + iPiF^ + 2jFB^ - H,pi - 3H,p2F^ - H^p^B^ = , (123) 

2aiA + UPiF^A + 27^5^ _ QH;,p2AF = -(1 + m)AD^ , (124) 

UpiFA"^ - 2m'yFB'^ - 3H^p2A'^ + mH^p^B'^ = , (125) 

2(3xA^ - m-yB"^ = mD'^ , (126) 

2a2B + {2/m)jA^B + 27^^^ - 2H^p^BF = (2 - m)BD'^ , (127) 

2mp2B^ - 7^2 = -mD^ , (128) 

AjFAB - 2HzpsAB = . (129) 

In this case also it turns out that the solution exists only if either 2(3i > 7 > 2/^2, 7^ > 4/3i/32 or 
2/^2 > 7 > 2/31,72 < APiP2- We find that 

2_ mD^[j-2P2] 2_ ^'[2/3i-7] 

[7'-4/3i/32] ' [7' -4/31/32]' ^ ^ 



where 



i2_o„d2 , 2pi7 P2P3H'^ iP_P3Hz 



-{l + m)D' = 2jB' + -^-^^^^, F = ^, (131) 

P3 27 27 



and the three constraints are 



^ P27 P2P3-f^2 , Pll 

Pi = 7^ ' "1 = o + . 

2p3 27 P3 

(2 - m)D'^ = 2a2 + {2/m)jA'^ - . (132) 



27 

Special case of Hz = 0: 
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With Hz = 0, the solution is again given by Eq. H122|) but with F = = G where A and B are again 
as given by Eq. (|13U|) and furthermore, ai,a2 turn out to be negative, i.e. 

ai = - \ ' iB^ , 02 = — - 2/32^2 . (133) 

m=l: 

In the special case of m = 1 and G = 0, -F 7^ 0, both solutions VI and VII reduce to a bright-dark type 
of solution as given by Eq. (|12U|1 . i.e. 

(/> = F + Atanh[D(x + xo)] , -0 = ^sech[L'(x + xq)] , (134) 
with A, B and D given by 

^ -[7^-W2]' -[7^-W2]' ''--^''-77 + ^^' ^'''^ 
while the other relations remain unchanged and are again given by Eqs. (|13H) and (|132jl . 
Special case of 7^ = 4/3i/32 

One can show that the solution l|122j) exists even in case 7^ = 4/3i/32. It turns out such a solution exists 
only if 

2/3i = 2/32 = 7, (136) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

A^ - mB^ = . (137) 

7 



Further, one has 



P3H! Pi7 p\Hl , , 

27 p3 2 27 



Energy: Corresponding to the mixed "kink-pulse lattice" solution [Eq. (|122j) ] the energy and the constant 
C are given by 

E = —( [(2m - 1)B^ + (1 + m)A'^]E(m) - (1 - m)(A'^ - B'^)K(m)] , 
3m \ J 

(139) 

C = -^J^D^ - F^[ai + SPiF^ - 2HzP2F] + [-7^^ + 02 + P2B^]B'^. (140) 
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For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbohc (bright-dark) sohton solution [Eq. ()12Ujl ] plus the interaction energy. We find 



E — Esoliton + Eint — D 



(141) 



.3 6 — m, 
Note that this solution also exists only when either 2(3i > 7 > 2/^2 and 7^ > 4/3i/32 or 2/^2 > 7 > 2/3i and 
7^ < 4/?i/32- The interaction energy vanishes for m = 1. 

2.8 Solution VIII 

Finally, we discuss two periodic solutions both of which at m = 1 reduce to a dark-bright type of solution. 
The first such (pulse-like in and kink-like in ip) solution is given by 

(j) = F + Acn[D{x + xo),m] , i/j = G + Bsn[D{x + xo),m] , (142) 

provided G = and the following seven coupled equations are satisfied 

2aiF + 4(3iF^ + 27^^^ _ ^^^^ _ 3H,p2F^ - H^p^B^ = , (143) 

2aiA + UpiF'^A + - QH^p2AF = {2m - 1)AD^ , (144) 

12/3iF^2 _ 2jFB^ - m,p2A^ + H.p^B^ = , (145) 

2piA^ - -fB^ = -mD^ , (146) 

2a25 + 2jA'^B + 2-fBF'^ - 2H^p^BF = -(1 + m)BD'^ , (147) 

2/?2-B2 - 7^2 ^ ^^2 ^ ^-^4g^) 

4-fFAB - 2H^p3AB = . (149) 

In this case it turns out that the solution exists only if either 2/?i > 7 > 2/^2, 7^ < 4/3i/32 or 2/?2 > 7 > 



2/3i, 72 > 4/?i/32. We find that 



[4/31/32-7'] ' [4/?i/?2-72] ' 

where 



(2m - 1)1)2 = 27^2 + ^ - , (151) 

P3 27 27 
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and the three constraints are 

Pi = T, ) «1 = 7, \ ■> 

-(l + m)D2 = 2a2 + 27A'-^^. (152) 

27 

Special case of Hz = 0: 

When = 0, the solution is again given by Eq. (|142jl but with F = = G where A and B are again 
as given by Eq. (|15(H) and furthermore, ai, a2 turn out to be negative, i.e. 

a^ = -^- 2i3,A^ , a2 = -il±|!^ _ ^A^ . (153) 

It is worth pointing out that with Hz = 0, the field equations are completely symmetric between the two 
fields (p and ip and hence the solutions VI and VIII are identical in the limit Hz = 0. 
Special case of 7^ = A[3i[32 

One can show that the solution H142() exists even in case 7^ = 4/3i/32. It turns out such a solution exists 
only if 

2/?i = 2/?2 = 7, (154) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

i?2 - ^2 = . (155) 

7 

Further, one has 

P2=P3, ai = ^ + — , «i-a2 = ^ + ^. 156 
27 /03 2 27 

Energy: Corresponding to the mixed "kink-pulse lattice" solution [Eq. (|142jl ] the energy and the constant 
C are given by 

E - Vmin = ^ {[{2m - l)A^ + (1 + m)B^]E{m) - (1 - m){B^ - A^)K{m)^ , (157) 

C = -^(1 - m)A^D'^ - F'^[ai + S^F^ - 2HzP2F] + [--/F'^ + 02 + P2B^]B'^. (158) 

For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-bright) soliton solution 

(j) = F + ^sech[L>(x + xq)] , ^p = B tanh[i:»(x + xq)] , (159) 
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plus the interaction energy. We find 



E — Egoliton + Eint — D 



1{2B^ + A^) + il-^(2i?2 - 5^2) + (1 - m)A^ In ( 

3 6 \\/l — mJ 



(160) 



Note that this solution also exists only when either 2(5i > 7 > 2/^2 and 7^ < 4/3i/32 or 2/?2 > 7 > 2/3i and 
7^ > 4/3i/?2. Again, the interaction energy vanishes at m = 1. It is amusing to note that the energy of the 
solution VIII is easily obtained from that of solution VII by simply interchanging A and B. 

2.9 Solution IX 

Another solution which at m = 1 reduces to a dark-bright type of soliton solution is given by 

<!) = F + Adn[D{x + XQ),m], = G + Bsn[D{x + XQ),m] , (161) 
provided G = and the following seven coupled equations are satisfied 

2aiF + 4/3iF3 ^ {2/m)^FB'^ - H,pi - 3H^p2F^ - ^^^^^ = , (162) 

2aiA + UPiF^A + {2/m)'yAB'^ - QH^p2AF = (2 - m)AD'^ , (163) 

UPiFA^ - (2/mhFB'^ - 3H,p2A'^ + = , (164) 

m 

2mPiA^ - 7^2 = -mD'^ , (165) 

2a25 + 2-^A^B + 2'yBF^ - 2H^pzBF = -(1 + m)BD'^ , (166) 

2/32^2 - mjA"^ = mD^ , (167) 

AjFAB - 2H^pzAB = . (168) 

In this case it turns out that the solution exists only if either 2/3i > 7 > 2/32, 7^ < 4/?i/32 or 2/^2 > 7 > 



2/3i, Y > 4/3i/32- We find that 



2_ iJ^[7-2/32] 2_ mD^Wi-l] 

[4/3i/32-7^]' [4/3i/32-7^] ' ^ ^ 



where 



(2 - m)D^ = (2/m)75^ + ^ - ^^^^ , = ^ , (170) 

P3 27 27 
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and the three constraints are 

-{l + m)D'^ = 2a2 + 2^A^ . (171) 

27 

Special case of = 0: 

When = 0, the solution is again given by Eq. (|161j) but with F = = G where A and B are again 
as given by Eq. (|169|) and furthermore, 01,02 turn out to be negative, i.e. 

mL>2 ^ „ (l+m)Z)2 „ ^ ^ 

ai = ^ - WiA^ , «2 = -^^y^ 1^ ■ (172) 

It is worth pointing out that with //^ = 0, the field equations are completely symmetric between the two 

fields (j) and ip and hence the solutions VII and IX are identical in the limit Hz = 0. 

m=l: 

In the special case of m = 1 and G = 0, F 7^ 0, both solutions VIII and IX reduce to a dark-bright 
type of solution H159|l . i.e. 

(j) = F + As<ic\i[D{x + xq)] , V = ^tanh[Z)(x + xo)] , (173) 
with A, B and D given by 

while the other relations remain unchanged and are given by Eqs. (|15H) and H152|l . 
Special case of 7^ = 4/3i/?2 

One can show that the solution H161() exists even in case 7^ = 4/3i/32. It turns out such a solution exists 
only if 

2/3i = 2/32=7, (175) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

- mA^ = . (176) 

7 
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Further, one has 



P2 = P3, ai = ^^ + ^, ai-a2 = — + ^^. 177 
27 p3 2 27 



Energy: Corresponding to the mixed "pulse-kink lattice" solution [Eq. H161jl ] the energy and the constant 
C are given by 

E - Vmin = ^ (^[(2 - m)mA^ + (1 + m)B'^]E{m) - (1 - m){B'^ + 2A^)K{m)^ , (178) 
1 

C=-{1- m)A^D'^ - F^[ai + SPiF^ - 2H,p2F] + [-m^F^ + maa + P2B^]^. (179) 

For m near one, the energy of this periodic solution can be rewritten as the energy of the corresponding 
hyperbolic (dark-bright) soliton solution [Eq. IjlSBj) ] plus the interaction energy. We find 



E — Esoliton + Eint — D 



6 6 \^/\ — mJ 



(180) 



Note that this solution also exists only when either 2/3i > 7 > 2/^2 and 7^ < 4/3i/32 or 2/^2 > 7 > 2/9i 
and 7^ > Aj3ij32- The interaction energy vanishes at m = 1. It is amusing to note that the energy of the 
solution IX is easily obtained from that of solution VI by simply interchanging A and B. 

Summarizing, we have obtained nine periodic solutions in terms of Jacobi elliptic functions, in the case 
of coupled (/>^ field theory with biquadratic coupling and an external magnetic field. This was possible 
because the magnetic field interaction is not symmetric between the two fields (j) and "0- In the special 
case when the modulus parameter m of the Jacobi elliptic function is one, these nine solutions reduce 
to four different soliton solutions valid on the full line and expressed in terms of hyperbolic functions. 
Note, however, that in case the external field Hz = 0, instead of nine, we only obtain six distinct periodic 
solutions (of which three are previously known [20) which in the limit m = 1, give three distinct soliton 
solutions. 

It is worth emphasizing the restrictions on the various parameters in the case of the nine solutions. 
For example, in the case of sn — sn solution (with G = 0), 2/?i > 7, 2/^2 > 7. Further, in the special case 
of Hz = = F = G, one can show that qi < 0, a2 < 0. On the other hand, in the case of cn — cn, 
dn — dn, cn — dn and dn — cn solutions (with G = 0), 7 < and further 7^ > 4/?i/32- In the special case 
oi Hz = = F = G, m addition one finds that (i) for cn — cn case ai, 02 > (<) provided m > (<) 1/2 
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(ii) «! > 0, a2 > for dn — dn solution (iii) for dn — cn solution, a2 > (iv) ai < for cn — dn solution. 
Instead, for sn — cn as well as sn — dn solutions, either 2/3i > 7 > 2/32 and 7^ > 4/?i/32 or 2p2 > 7 > Wi 
and 7^ < 4/32/34- Finally, in the cn — sn and dn — sn cases, either 2/3i > 7 > 2/^2 and 7^ < 4/?i/32 or 
2/^2 ^ 7 ^ 2/?! and 7^^ > Aj3ij32- If in addition = = F = G, then for all the four solutions (i.e. sn — cn, 
cn — sn, sn — dn, dn — sn), qi < and 02 < 0. 



3 Solutions with bilinear coupling 

Several years ago, a coupled (/)^ model was considered in the context of a surface phase transition with 
hydration forces which is similar to the one considered in the last section except that there was no 
external magnetic field Hz and instead of a biquadratic coupling there was a bilinear coupling between the 
two fields. The purpose of this section is to obtain a bright-bright soliton solution of that model. 
The potential (i.e. free energy) is given by 

V = ai(t? + /?i</.^ + a2^^ + /32^^ + 5i{<P- ^pf + ^2(0 + V')^ , (181) 

with (5i, 82 being the coupling parameters between the two fields. The (static) equations of motion which 
follow from here are 

(Pd> 

= 2ai0 + 4/3103 + 25i{^ - 0) + 2<52(0 + ^) , (182) 
^ = 2a2^p + 4/?2V'^ - 25i((/) - V') + 262{(l) + ^) . (183) 

It is not difficult to show that this pair of coupled field equations admits the "kink-kink" type periodic 
solution 

(j) = Asn[D{x + XQ),m] , ip = Bsn[D{x + XQ),m] , (184) 

provided 

mD^ = 2/3iA^ = 2/32^^ , 



a2 



(1 + m)/3i 



01 + 81 + 62 + {52 - Si)J — 



(185) 
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and furthermore the parameters ai,a2, f3i, f32,Si,52 satisfy the constraint 

a2 — ai = {Si — 62) 



IP2_ (h 
Pi V/32 



(186) 



Since > 0, the relation (|185j) gives us a strong constraint on some of the parameters. The energy E 
and the constant C corresponding to the periodic solution [Eq. (|184jl ] are 

12 , Ty2\ 



3m 



[(1 + m)E{m) - (1 - m)K{m)] , 



1 



C = --{A'^ + B^)D^ . 



In the limit m = 1, this solution reduces to the bright-bright soliton solution 



(187) 



Atanh[D{x + xo),m] , ip = B tanh[D{x + xq)] 



(188) 



provided 



A' 



-1 



ai + 81 + 82 + {82 - 



(189) 



while the relation H186|) remains unchanged. For m near one, the energy of the periodic solution can be 
rewritten as the energy of the corresponding hyperbolic (bright-bright) soliton solution [Eq. H188(l ] plus 
the interaction energy. We find 



E = Ekink + Eint = {A'^ + B'^)D 



4 (1 - m) 

- + 

3 3 



(190) 



The interaction energy vanishes at m = 1. In view of the requirement Pi, P2 > arising from stability, 
we are unable to find any other solution to this coupled set of equations with a bilinear coupling. 



4 Solutions of discrete coupled 0*-type equations with biquadratic cou- 
pling 

Discrete coupled (j)^ models arise in the context of structural transitions on a lattice, collective proton 
dynamics in ice [23]) etc. The purpose of this section is to give an exhaustive list of solutions to the 
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discrete coupled -type equations with biquadratic coupling (but in the absence of an external magnetic 
field Hz). In the next section, we shall obtain a solution of the discrete coupled c^^-type equations with 
bilinear coupling. 

We start from the coupled static field equations and ©• The discrete analog of these field equations, 
for Hz = has the form 

^(0„+i + - 2</.„) - 2qi(/>„ - 2[2Ml + 7V'^]0„ = , (191) 

^(V-n+l + V-n-l - 2V'„) - 2a2i^n " 2[2/32^^ + T'/'nlV'n = . (192) 

We are unable to find any solution to this coupled set of field equations. However, as in the Ablowitz-Ladik 
discretization of the discrete nonlinear Schrodinger equation 24 , if we replace (j)n and ipn in the last term 
in Eqs. ()191|) and (|192() by their average, then we can find exact solutions to this coupled system. In 
particular, instead of Eqs. (|191() and H192() . we consider the discretized equations 

^{(f>n+l + (f>n-l - 2(f>n) " 2ai0„ - [2Pi(t>l + li^lUn+l + 4>n-l] = , (193) 
■^{ipn+l + i'n-l - 2ll>n) - 2a2'4^n - [2MI + 7</'n][V'n+l + ^n-l] = . (194) 

Note that single solitons and their stability in coupled Ablowitz-Ladik chains have been studied previously 
j25j . We now show that this modified set of coupled discrete equations has six different periodic solutions 
which in the limit m = 1 reduce to the bright-bright, bright-dark and dark-dark soliton solutions. In all the 
solutions, we shall see that the static kink can be placed anywhere with respect to the lattice. Hence we 
suspect that in all these cases, there may be an absence of the Peierls-Nabarro barrier [261 1271 05] . which 
is the energy cost associated with moving a localized solution such as a soliton by a half lattice constant 
on a discrete lattice. It would be nice if one can demonstrate this explicitly. 

4.1 Solution I 

It is easy to show that the field Eqs. H193|) and H194|) admit the kink-kink type solution 

(fin = -A- sn[hD{n + XQj^m] , Vn = sn[/iD(n + xq), , (195) 
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provided 



2 m(2/32 — 7)sn^(/i-D, m) 2 m(2/?i — 7)311^ (/iD, m) 
/i2(4/3i/32-72) ' " /i2(4^^^2-7 



2^ 



Oi = "2 = ~J^i^ ~ cn(/iD, m)dn(ft.L', m)], (196) 

where /i is the lattice spacing. Thus, note that as in the continuum case, this solution exists provided 
2/?! > 7, 2/^2 > 7, ai < 0, 02 < 0. It is interesting to note that the solutions to both the discrete and the 
continuum model exist under the same set of conditions. 

Continuum Limit: It is instructive to consider the continuum limit h ^ and show that the above 
solution smoothly goes over to the corresponding continuum solution. In particular, on using the fact that 
as /i ^ 

sn{hD, m) ^hD, cn^{hD, m) ^ I - h^D"^ , dn^{hD, m) I - mh'^D^ , (197) 

it readily follows that the above solution indeed reduces to the corresponding continuum solution [Eq. (7)] 

obtained in Sec. II (when = F = G = 0), i.e. 

2 m(2/32-7)^' „2 m{2f3,-j)D^ {l + m)D^ 

(4/3i/32 - 7 ) (4/3i/32 - 7 ) 2 

In fact we shall see that all the six solutions of this coupled discrete model smoothly go over to the 
corresponding continuum solutions obtained in Sec. II in the limit h ^ 0. 

In the limit m = 1, the periodic solution H195|) reduces to the bright-bright soliton solution 

(j)n = A tanh[hD{n + xq)] , t/^n = B tanli[hD{n + xq)] . (199) 

Special case of 7^ = 4/?i/?2 

One can show that the solution (|195|) exists even in case 7^ = 4/3i/32. It turns out that such a solution 
exists only if 

2/3i = 2/32 = 7 , (200) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

A2 + 52 = ^^^^4?^- (201) 

In the continuum limit /i — > 0, as expected, this reduces to the constraint equation obtained in Sec. 
II. 
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4.2 Solution II 

It is easily shown that a kink-pulse type solution 

(l)n = A sn[hD{n + xo),m] , ipn = B cn[hD{n + xo),m] , (202) 

is an exact solution to the field Eqs. H193() and H194() provided 

2 ^ m{-f -2f32)sn'^{hD,m) . ^ m(2/3i - j)sn'^ {hP , m) 

Furthermore, 

cn{hD,m)[2m(3i{j - 2(52WihD,m) - (7^ - 4/?i/?2)] 
/i2 + /i2(72 -4/3i/52)dn(/iZ),m) ' ^ ' 

cn(^D,m)[m7(7 - 2(32W{hD,m) - (72 - 4/3i/?2)] 
/i2+ /i2(72_4/3i/32)dn2(/iD,m) ' ^ ' 

Again, as in the continuum case either 2/?i > 7 > 2/?2, 7^ > 4/?i/32 or 2(52 > 7 > 2/3i, 4/3i/32 > 7^- It is, 

however, not clear here if ai and 02 have a definite sign. However, in the limit m = 1, as in the continuum 

case, one finds that ai < 0, Q2 < 0. 

It is readily checked that in the continuum limit this solution smoothly goes over to the corresponding 
continuum solution, Eqs. (104) and H112|) . Furthermore, the corresponding bright-dark solution is easily 
obtained in the limit m = 1. 
Special case of 72 = 4/3i/?2 

One can show that the solution (|2()2j) exists even in case 7^ = 4/?i/?2. It turns out such a solution exists 
only if 

2/3i = 2/32 = 7 , (206) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

- B^dn^ihD, m) = ,\ ' ^ . (207) 

n^7 

In the continuum limit ^ — > 0, as expected, this reduces to the constraint equation (|117|) obtained in Sec. 
II. 
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4.3 Solution III 

Yet another kink-pulse type solution is given by 

(j)n = A sn[/iZ)(n + xq), m] , ipn = B dn[hD{n + xo),m\ , (208) 



provided 



,2 m(7- 2/32)sn2(/iD,m) . (2^1 - j)sn^ihD,m) 



Furthermore, 



_ 1 dn(/iZJ,m)[2/?i(7-2/32)sn^(/^I^,m)-(7^-4/3i/32)] 

/i2 + /i2(72 -4/3i/32)cn(/iD,m) ' ^ ^ 

_ 1 dn(/iA^)[7(7-2/?2)sn^(/^Zj,m)-(7^-4A/32)] . . 

/i2 ^ /i2(72 -4/3i/32)cn2(/iD,m) " ^ ' 

Special case of 7^ = 4/3i/?2 

One can show that the solution (208) exists even in case 7^ = 4/3i/52. It turns out such a solution exists 
only if 

2/?i = 2/?2 = 7 , (212) 
and that in this case one cannot determine A, B. However, they must satisfy the constraint 

- mB^cn^(hD, m) = ,\ ' ^ . (213) 

/i^7 

In the continuum limit /i — > 0, as expected, this reduces to the constraint equation (|137|) obtained in Sec. 
II. 

Note that in the m = 1 limit the solutions II and III reduce to the same bright-dark soliton solution. 
In addition, as in the continuum case, this solution exists if either 2/?i > 7 > 2/?2, 7^ > 4/3i/32 or 
202 > 7 > 2/?!, 4/3i/?2 > 7^. It is, however, not clear here if qi and 02 have a definite sign. However, in 
the limit m = 1, as in the continuum case, one finds that qi < and a2 < 0. 

4.4 Solution IV 

Finally, we present three periodic solutions, all of which in the limit m = 1 reduce to the dark-dark soliton 
solution. One of the pulse-pulse type periodic solution is given by 

4>n = A cii[hD{n + xo), m] , tpn = B cn[hD{n + xo),m] , (214) 
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provided as in the continuum case, 7 < and 7 > 4/3i/32. We find 

,2 2m(/?2 + |7l)sn2(/iZ),m) „ 2m(/?i + |7|)sn2(/iZ), m) 

A = — — — ^ rr- , i) — 



/i2(^2 _ 4/3i4/?2)dn2(/iD, m) ' /i2(72 - 4/3i4/32)dn2(/iD, m) 

1 



^ cn{hD,m) 



(215) 



dn^ (/iL»,m) 

Using Eq. ()197() it is easily shown that as in the continuum case, ai, Q2 > (<) provided m > (<) 1/2. 
This solution is equivalent to the continuum solution, Eqs. (33) and (42). 

4.5 Solution V 

Another pulse-pulse type solution is given by 

(pn = ^ dn[hD{n + xo),m] , tpn = B dn[hD{n + xq), m] , (216) 

provided as in the continuum case, 7 < and 7^ > 4/3i/32. We find 

,2 2{p2 + h\)sn^hD,m) . 2{(3, + \-f\)sn\hD,m) 



/i2(^2 _ 4/3i4/32)cn2(/iD, m) ' h?{-i'^ - 4/3i4/32)cn2(/iZ), m) ' 

1 



ai = "2 = -j^ 



^ dn(/iL',m) 



(217) 



cn2(/iL), m) 

Using Eq. H197|) it is easily shown that as in the continuum case, ai, 02 > 0. This solution is equivalent 
to the continuum solution, Eqs. (33) and (62). 

4.6 Solution VI 

Yet another pulse-pulse type solution is 

(pn = A dn[hD{n + xo),m] , tpn = B cn[hD{n + xo),m] , (218) 

provided as in the continuum case, 7 < and 7^ > 4/3i/32- We find 

^2 ^ i2(32 + h\)sn\hD,m) ^2 ^ m(2/?i + |7|)sn2(/,ZJ, m) 

/i2(72 - 4/3i^2)cn2(/iZ),m) ' /i2(^2 _ 4.Pip2)dn\hD,m) ' ^ ' 



Furthermore, 



^ 1 , mm + |7l)dn"(/^^,^) - 7(2/?i + |7l)cn2(/iZj,m))] 

/i2(72 _4/3i/32)dn(/iZ),m)cn2(/iZ),m) ' ^ ^ 
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_^ ^ 1 W2{2l3i + \-i\W{hD,m)-\^m+2(52)dn\hD,m)] 

h? h^{-f^ -APiP2)dn^(.hD,m)cn{hD,m) ' ^ ^ 

It is not clear if in general ai and 02 have a definite sign. However, it is easily checked that at m = 1, ai, 

a2 > 0. This solution is equivalent to the continuum solution, Eqs. (72) and (80). 

Note that the last three (i.e. IV,V,VI) solutions reduce to the (same) dark-dark soliton solution in the 

limit m = 1. Since we do not know the Hamiltonian corresponding to Eqs. (193) and (194), we are unable 

to find the energy and soliton interaction explicitly for any of the above discrete solutions. Similarly, for 

Hz / we have not succeeded in finding exact solutions. 

5 Solutions of discrete coupled 0^-type equations with bilinear coupling 

Coupled lattice chains, with a bilinear coupling, undergoing a second order structural phase transition can 
represent this case. We start from the coupled static field Eqs. (|182() and H183() . The discrete analog of 
these field equations has the form 

^{(Pn+l + <Pn-l - 2<Pn) " 2[qi + <5i + Js]^™ + 2[6i - dsJV'n " ^l3i<pl = , (222) 

^(V'n+l + V'n-l - 2Vn) " 2[a2 + 5l + <52]^n + 2[<5i - 62]<l)n " 4/32^^ = . (223) 

We are unable to find any solution to this coupled set of field equations. However, as in the Ablowitz-Ladik 
discretization of the discrete nonlinear Schrodinger equation [24., if we replace 0„ and ipn in the last term 
in Eqs. (|222|) and ()223() by their average, then we can find exact solutions to this coupled system. In 
particular, instead of Eqs. ()222() and H223() . we consider the discretized equations 

^((An+l + (pn-l - 2(Pn) - 2[ai + 5i + J2]</'n + 2[5i - 52]'4>n - 2f3i<i)l[^n+l + 0n-l] = , (224) 
^(V'n+l + V'n-l - 2ll)n) - 2[a2 + 6i + ^2]^'^ + 2[5i - 52](t)n - 2/^2 [V'n+1 + V'n-l] = . (225) 



It is easy to show that the field Eqs. (|224j) and (|225|) admit the kink-kink type solution 

= j4 sii[hD{n + xo), m] , ipn = B sn[hD{n + xo),m] , (226) 
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provided 



msn'^{hD, m) 



2/3i[l - cn{hD,m)<ln{hD,m)] 



—•msT?{hD, m) 



ai + 5i + 62 + {82 - Si) 




(227) 



and furthermore the parameters ai,a2, Pi, /32,Si, 62 



satisfy the constraint 



"2 — ai = (<5i — h) 




(228) 



As expected, in the continuum hmit of /i ^ 0, this solution smoothly goes over to the corresponding 
continuum solution, Eqs. (184), (185) and (186), obtained in Sec. III. Since we do not know the Hamiltonian 
corresponding to Eqs. (224) and (225), we are unable to find the energy and soliton interation explcitly 
for this discrete solution. Similarly, for 7^ we have not succeeded in finding an exact solution. 

6 Conclusion 

We have systematically provided an exhaustive set of exact periodic domain wall solutions for a coupled 0^ 
model with and without an external field, and for both bilinear and biquadratic couplings. Only a bright- 
bright solution could be obtained for the bilinear case. For both the biquadratic and bilinear couplings the 
corresponding discrete case was also considered-with an Ablowitz-Ladik like modification of the coupled 
discrete equations- and we obtained several exact solutions. For the solutions of the discrete model, the 
calculation of the Peierls-Nabarro barrier |261 1271 128j and soliton scattering |291 l3Uj remain topics of further 
study. Similarly, scattering of solitons in the coupled (/)^ continuum and discrete models with either the 
biquadratic or bilinear coupling is an interesting open issue. To this end, the static solutions presented 
here need to be boosted with a certain velocity. 

It would be instructive to explore whether the nine different solutions reported in Sec. 2 (or the six 
solutions in Sec. 4) are completely disjoint or if there are any possible bifurcations linking them via, for 
instance, analytical continuation. We have not tried to carry out an explicit stability analysis of various 
periodic solutions. However, the energy calculations and interaction energy between solitons (for m ~ 1) 
in the case of both the biquadratic and bilinear couplings could provide useful insight in this direction. 
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Our results are relevant for spin configurations, domain walls and magnetic phase transitions in mul- 
tiferroic materials [SlHjj periodic domain walls are yet to be observed in the hexagonal multiferroics 
Similarly, our solutions are important for understanding structural phase transitions in ferroelectrics E] 
and elastic materials biophysics problems such as multilamellar lipid systems as well as field theo- 
retic contexts j^lEl- These ideas and exact solutions can be generalized to other coupled models such as 
(p^ (for first order phase transitions) and will be discussed elsewhere j31j . 

7 Acknowledgment 

A.K. acknowledges the hospitality of the Center for Nonlinear studies at LANL. This work was supported 
in part by the U.S. Department of Energy. 

References 

[1] H. Schmid, Ferroelectrics 162, 317 (1994); N. A. Spaldin and M. Fiebig, Science 309, 391 (2005). 

[2] T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima, and Y. Tokura, Nature 426, 55 (2003). 

[3] M. Fiebig, Th. Lottermoser, and R. V. Pisarev, J. Appl. Phys. 93, 8194 (2003). 

[4] S. H. Curnoe and I. Munawar, Physica B 378-380, 554 (2006). 

[5] S. Aubry and R. Pick, Ferroelectrics 8, 471 (1973). 

[6] T. Abel and R. Siems, Ferroelectrics 153, 177 (1994). 

[7] A. A. Kornyshev, D. A. Kossakowski, and S. Leikin, J. Chem. Phys. 97, 6809 (1992). 

[8] A. N. Das and B. Ghosh, J. Phys. C 16, 1803 (1983). 

[9] R. Rajaraman, Phys. Rev. Lett. 42, 200 (1979). 

[10] C. Matsuoka and K. Nozaki, Phys. Lett. A 185, 310 (1994). 

[11] C. S. Lai, Canad. J. Phys. 58, 443 (1980). 

34 



[12] N. N. Rao, J. Phys. A 22, 4813 (1989). 

[13] X. Y. Wang, B. C. Xu, and P. L. Taylor, Phys. Lett. A 173, 30 (1993). 

[14] X. W. Huang, J. H. Han, K. Y. Qian, and W. Qian, Phys. Lett. A 182, 300 (1993). 

[15] D. Bazeia, M. J. Dos Santos, and R. F. Ribeiro, Phys. Lett. A 208, 84 (1995). 

[16] Y. J. Zhu and S. Y. Lou, Commun. Theor. Phys. 30, 147 (1998). 

[17] S. Y. Lou, J. Phys. A 32, 4521 (1999). 

[18] D. B. Cao, Phys. Lett. A 296, 27 (2002). 

[19] D. S. Li and H. Q. Zhang, Acta Phys. Sinica 52, 2373 (2003); ibid. 2379 (2003). 

[20] C. P. Liu, Chaos Soht. Fract. 20, 619 (2004); Commun. Theor. Phys. 43, 13 (2005). 

[21] S. K. Liu, Z. T. Fu, S. D. Liu, and Z. G. Wang, Phys. Lett. A 323, 415 (2004). 

[22] Handbook of Mathematical Panctions, eds. M. Abramowitz and L A. Stegun (Dover, New York, 1972). 

[23] A. V. Zolotaryuk, A. V. Savin, and E. N. Economou, Phys. Rev. B 57, 234 (1998). 

[24] M. J. Ablowitz and J. F. Ladik, J. Math. Phys. 16, 598 (1975); ibid. 17, 1011 (1976). 

[25] B. A. Malomed and J. Yang, Phys. Lett. A 302, 163 (2002). 

[26] F. R. N. Nabarro, Theory of Crystal Dislocations, (Dover, New York, 1987). 

[27] O. M. Braun and Yu. S. Kivshar, Phys. Rev. B 43, 1060 (1991). 

[28] P. G. Kevrekidis, Physica D 183, 68 (2003). 

[29] A. Khare, K. 0. Rasmussen, M. R. Samuelsen, and A. Saxena, J. Phys. A: Math. Gen. 38, 807 (2005); 
A. Khare, K. 0. Rasmussen, M. Salerno, M. R. Samuelsen, and A. Saxena, Phys. Rev. E 74, (2006). 

[30] F. Cooper, A. Khare, B. Mihaila, and A. Saxena, Phys. Rev. E 72, 036605 (2005). 

[31] A. Khare and A. Saxena, to be published. 

35 



